Medium-mass nuclei from chiral nucleon-nucleon interactions 
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We compute the binding energies, radii, and densities for selected medium-mass nuclei within 
coupled-cluster theory and employ the "bare" chiral nucleon-nucleon interaction at order N 3 LO. We 
find rather well-converged results in model spaces consisting of 15 oscillator shells, and the doubly 
magic nuclei 40 Ca, 48 Ca, and the exotic 48 Ni are underbound by about 1 MeV per nucleon within 
the CCSD approximation. The binding-energy difference between the mirror nuclei 48 Ca and 48 Ni is 
close to theoretical mass table evaluations. Our computation of the one-body density matrices and 
the corresponding natural orbitals and occupation numbers provides a first step to a microscopic 
foundation of the nuclear shell model. 
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Introduction. Ab-initio nuclear structure calculations 
have made great progress in the past decade. Light nu- 
clei up to carbon or so can now be described in terms of 
their nucleonic degrees of freedom and realistic nucleon- 
nucleon (NN) forces (i.e. those that include pion ex- 
change and fit the NN phase shifts up to 350 MeV lab 
energy with a x 2 ~ l_P er datum) augmented by a three- 
nucleon force (3NF) 0,0, 0. One of the major advances 
is due to the systematic construction of nuclear forces 
within chiral effective field theory (EFT) 0, @ . In this 
EFT, unknown short-ranged physics of the nuclear force 
is systematically parametrized in terms of contact terms 
and their low-energy constants, while the long-range part 
of the interaction stems from pion exchange. One of the 
hallmarks of this approach is the "power counting", i. 
e. an expansion of the nuclear Lagrangian in terms of 
the momentum ratio Q/A. Here, Q denotes the typical 
momentum scale at which the nucleus is probed, while A 
denotes the high-momentum cutoff scale that limits the 
applicability of the EFT. Within this approach, three- 
nucleon forces appear naturally at order (Q/A) 3 and 
four-nucleon forces appear at order (Q/A) 4 (EHS]. 

The chiral interactions have been probed in light sys- 
tems up to mass 13 
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Fujii et al. have 

employed chiral NN interactions for studies of 16 O [l3| 
within the unitary- model-operator approach (UOMA). 
Unfortunately, virtually nothing is known about chiral 
interactions in heavier nuclei. In particular, a study of 
their saturation properties is missing, and the contribu- 
tions of chiral NN interactions to nuclear binding and 
structure in medium-mass nuclei needs to be determined. 
It is the purpose of the present Letter to fill this gap. 

Ab initio methods began to explore medium-mass nu- 
clei only very recently. Gandolfi et al. 1A\ employed the 
auxiliary field diffusion Monte Carlo method for a com- 
putation of the binding energy of 40 Ca. However, this 
impressive calculation is not entirely realistic since the 
employed Argonne Vq potential lacks the spin-orbit in- 
teraction. Roth and Navratil [IH employed softer renor- 



malized NN interactions and computed the binding en- 
ergy of Ca within an importance truncated no-core 
shell model approach. However, this calculation was crit- 
icized [161.1171 for its convergence properties, the violation 
of Goldstone's linked cluster theorem and the correspond- 
ing lack of size extensivity. In this Letter, we employ the 
"bare" chiral NN interaction [3] and employ the size ex- 
tensive coupled-cluster method [II, [11, 0, [H HI H3 for 
the computation of various properties of the medium- 
mass nuclei 40 Ca, 48 Ca, and the exotic 48 Ni. The use 
of the "bare" NN interaction has the advantage that it 
avoids the introduction of additional many-body forces 
that are typically generated in secondary renormalization 
procedures of the two-body force. While our calculation 
includes the chiral NN interaction [?[ at order N 3 LO, it 
neglects the contributions of any 3NFs. 

This Letter is organized as follows. First, we briefly 
introduce spherical coupled-cluster theory. Second, we 
compute the binding energies of the nuclei 4 He, 16 0, 
40 Ca, 48 Ca, and 48 Ni with the "bare" chiral NN poten- 
tial. 

Spherical coup led- cluster theory. Coupled-cluster the- 
ory [3 H3 ED, HI 01 is based on the similarity 
transform 



H = e~ T He T 
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of the normal-ordered Hamiltonian H. Here, the Hamil- 
tonian is normal-ordered with respect to a product state 
\4>) which serves as a reference. Likewise, the particle- 
hole cluster operator 



T = T 1 +T 2 + . 



(2) 



is defined with respect to the reference state. The k- 
particle fc-hole (fcp-fch) cluster operator is 



T k = 
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Here and in the following, i, j, k, . . . label occupied single- 
particle orbitals, while a,b,c, . . . label unoccupied or- 
bitals of the reference state, i.e. it should have signif- 
icant overlap with the ground state. Throughout this 
work we will restrict ourselves to the CCSD approxima- 
tion f « f x +f 2 . The unknown amplitudes tf and in 
Eq. ([2]) are determined from the solution of the coupled- 
cluster equations 
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Here |#> 



is a lp-lh excitation of the reference 



state, and \<t>1j) is a similarly defined 2p-2h excited state. 
The CCSD equations ^ and ([5]) thus demand that the 
reference state \4>) is an eigenstate of the similarity trans- 
formed Hamiltonian (fTJ) in the space of all lp-lh and 2p- 
2h excited states. Once the CCSD equations are solved, 
the correlation energy of the ground state is computed 
as 



\H\ 



(6) 



Coupled-cluster theory fulfills Goldstone's linked clus- 
ter theorem and therefore yields size-extensive results. 
This is particularly important in applications to medium- 
mass nuclei. Within the CCSD approximation, the com- 
putational effort scales as n^n 4 where n a and n u denote 
the occupied and unoccupied orbitals of the reference 
state respectively. Thus, the computational effort 
is much smaller than within the configuration interac- 
tion for a given model space. This method has recently 
been employed in several ab initio nuclear structure cal- 
culations [2J, |25j, It is also ab e to compute 
lifetimes of unstable nuclei 2j| , to treat 3NFs [13] , and 
it meets benchmarks [l7| . 

For spherical reference states (i.e. nuclei with closed 
major shells or closed subshells), one can employ the 
spherical symmetry to further reduce the number of un- 
knowns (i.e. the number of cluster amplitudes). For such 
nuclei, the cluster operator ([2]) is a scalar under rotation, 
and depends only on reduced amplitudes. A naive esti- 
mate shows that a model space of n a + n u single-particle 
states consists of only (n G + n u ) 2 / 3 j-shells. Thus, the 
entire computational effort is approximately reduced by 
a power 2/3 within the spherical scheme compared to 
the m-scheme. We have derived and implemented the 
spherical scheme within the CCSD approximation. We 
tested that our m-scheme code and the spherical code 
give identical results for several cases. 

Results. The single-particle basis consists of wave 
functions of the spherical harmonic oscillator with the 
spacing hu>, the radial quantum number n, and angular 
momentum I, and we include single-particle states with 
In + I < N in our model space. The largest model space 
we consider (N = 14) consists of 15 oscillator shells. In 



such a large model space, configuration interaction be- 
comes impossible as the proton space alone consists of 
about 10 40 Slater determinants for 40 Ca. We first trans- 
form the Hamiltonian to the spherical Hartree-Fock ba- 
sis, and the CCSD equations are solved in this basis. 
Fully converged observables must be independent of the 
parameters N and fku of our single-particle basis. In 
practice, we cannot go to infinitely large spaces, and the 
dependence of our results on these parameters serve to 
gauge the convergence. 
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FIG. 1: (Color online) CCSD binding energy for 4 He from 
the chiral NN potential at order N 3 LO as a function of the 
oscillator spacing hui and the size of the model space. 

As a test case, Fig. [1] shows that the CCSD results for 
4 He are converged within a few keV with respect to in- 
creases in the size of the model space (denoted by N) and 
variation of the oscillator frequency. The triples correc- 
tions are not yet available within the spherical scheme, 
and we employ our m-scheme code for this purpose. The 
CCSD-T1 triples correction [31] in model spaces up to 
N = 7 yields another 1.3 MeV of binding. Thus, the 
CCSD-T1 results are very close to the virtually exact 
Faddeev-Yakubowski result E — —25.41 MeV quoted in 
Ref. for the same chiral NN interaction. The ex- 
perimental value is E = —28.3 MeV, and the additional 
binding is due to the missing 3NFs. 

The CCSD energies for 16 (see Fig. [2]) are converged 
within the order of about 100 keV and change by less than 
1 MeV over a considerable variation of the oscillator fre- 
quency. This result is in reasonably good agreement with 
the work by Fujii et al. who obtained -110 MeV as the 
binding energy from the UOMA 13]. Recall that both 
methods are approximations and based on similarity- 
transformed Hamiltonians. 

We turn to nuclei in the mass-40 region. The CCSD 
results for 40 Ca are shown in Fig. [3] Increasing the model 
space from TV = 13 to N = 14 yields an additional 
0.9 MeV, and the fiu-dependence is less than 2.2 MeV 
over the considered range of oscillator frequencies. Thus, 
the convergence with respect to the parameters of our 
model space is very satisfactory, and we are missing about 
10% of the experimental binding energy of -342MeV. 
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FIG. 2: (Color online) CCSD binding energy for 16 O from 
the chiral NN potential at order N 3 LO as a function of the 
oscillator spacing hw and the size of the model space. 
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FIG. 4: (Color online) CCSD binding energy for 48 Ca from 
the chiral NN potential at order N 3 LO as a function of the 
oscillator spacing hui and the size of the model space. 
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FIG. 3: (Color online) CCSD binding energy for 40 Ca from 
the chiral NN potential at order N 3 LO as a function of the 
oscillator spacing fku and the size of the model space. 



We also computed the binding energy of the mirror 
nuclei 48 Ca and 48 Ni. For 48 Ca the convergence of the 
results is satisfactory as shown in Fig. 2J and the con- 
vergence is very similar for 48 Ni. 48 Ni was discovered 
only recently [32j. It is believed to be a two-proton 
emitter, and its lifetime is very large compared to a 
typical nuclear time scale (i.e. the "orbital period" of 
a nucleon inside the nucleus). Thus, we can describe 
48 Ni in terms of a spherical Hartree-Fock basis based 
on the oscillator orbitals. Recall that the chiral inter- 
action includes charge symmetry-breaking and charge 
independence-breaking effects, and we also included the 
Coulomb interaction. The difference of our CCSD results 
for the mirror nuclei 48 Ca and 48 Ni is 1.38 MeV per nu- 
cleon and stems from these combined effects. Theoretical 
mass table evaluations 33;] suggest that the binding en- 
ergy of 48 Ni is 1.43 MeV per nucleon smaller than for 
48 Ca. Our results are in good agreement with this es- 
timate. The density of 48 Ca is shown in Fig. [5] The 
results still exhibit a dependence of the oscillator spac- 
ing huj, and the central density decreases with decreasing 
hw. This observable is less well converged than the en- 



ergy with respect to the size of the model space. The 
convergence is slow with respect to the maximum radial 
quantum number n employed in our model space, while 
the single-particle angular momentum I could be limited 
to I < 7. 




FIG. 5: (Color online) Densities for 48 Ca from the chiral NN 
potential at order N 3 LO for four different values of the oscil- 
lator spacing hw. 

Table U summarizes some of our results which are taken 
at fiuj = 28MeV in the largest model spaces. We com- 
puted the potential energy V via the Hellman-Feynman 
theorem. The fourth column shows the energy deviation 
AE = E — i?oxp from the experimental binding energy 
E exp for the considered nuclei. This difference is mainly 
due to the omitted 3NFs and the missing triples correc- 
tion. Note that 40 Ca is particularly tightly bound when 
compared to the other nuclei. The isotopes 16 0, 48 Ca, 
and 48 Ni all lack about AE/A ks 1.2 MeV of binding en- 
ergy when compared to experiment, while this difference 
is considerably smaller for 40 Ca. This result is somewhat 
surprising since 48 Ca is thought to be a better example 
of a doubly magic nucleus than 40 Ca. There seems to 
be a cancellation between triples corrections and contri- 
butions of 3NFs in 40 Ca. In other words, the isospin- 
dependence and/or mass-dependence of the 3NF is ex- 



4 



pected to be non-trivial. The charge radii are corrected 
according to Ref . [34| to account for the finite charge radii 
of the nuclcons. They are computed from the leading 
approximation of the center-of-mass corrected intrinsic 
density (25|. Note that the radii change about 0.1-0.25 
fm as the oscillator spacing hu> is varied in the range that 
is shown in the previous figures, and they decrease with 
increasing values of huj. 



Nucleus 


E/A 


V/A 


AE/A 


R 


Rcxp 


4 He 


-5.99 


-22.75 


1.08 


1.86 


1.64 


16 Q 


-6.72 


-30.69 


1.25 


2.71 


2.74 


40 Ca 


-7.72 


-36.40 


0.84 


3.24 


3.48 


48 Ca 


-7.40 


-37.97 


1.27 


3.22 


3.47 


48 Ni 


-6.02 


-36.04 


1.21 


3.50 





TABLE I: CCSD results for various nuclei from the chiral 
N 3 LO nucleon-nucleon potential. The binding energy per nu- 
cleon, and potential energy per nucleon are denoted as E/A 
and V/A, respectively. AE denotes the difference to the ex- 
perimental binding energy (difference to theoretical mass ta- 
ble evaluations for 48 Ni). R and i? C xp denote the computed 



and measured charge radius, 
and lengths in units of fm. 



Energies are in units of MeV, 



We also compute the one-body density matrices p pq — 
-Jo g ) of the ground states within the equation-of-motion 
CCSD [35] . The diagonalization of this matrix yields nat- 
ural orbitals and the corresponding occupations. These 
model-dependent quantities are, of course, not observ- 
ables but rather tied to the specific interaction we em- 
ployed. The dominant occupation probabilities are larger 
than 0.95, and this indicates that the considered nuclei 
are indeed doubly magic. This result is non-trivial. Note 
that the Hartree-Fock approximation does not even yield 
bound nuclei. Yet the CCSD correlations imprinted onto 
the Hartree-Fock state yield a rather simple state. To 
our knowledge, this is the first time the phenomenolog- 
ical shell-model picture of independent nucleon motion 
arises within an ab-initio approach. 

Summary. We have studied the saturation properties 
of chiral NN interactions at the order N 3 LO in medium- 
mass nuclei within the CCSD approximation of coupled- 
cluster theory. Our results exhibit a very satisfactory 
convergence with respect to the size of the model space 
and are only weakly dependent on the oscillator param- 
eter. We find that the "bare" chiral NN potential un- 
derbinds nuclei by about 1 MeV per nucleon. The com- 
parison of 40 Ca with 48 Ca and 48 Ni hints at an isospin 
dependence of the 3NF in medium-mass nuclei. Within 
the CCSD approximation, the proton-rich nucleus 48 Ni 
is less tightly bound by f .38 MeV per nucleon than its 
mirror nucleus 48 Ca, and this result is in good agreement 
with theoretical mass table evaluations. These calcula- 
tions pave the way to probing chiral interactions in even 
heavier nuclei and link the phenomenological shell model 



to ab-initio calculations. 
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